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1.  INTRODUCTION 

The  flow  in  centrifugal  pumps  is  exceedingly  complex,  involving  curvature,  system  rotation,  separation, 
turbulence,  and  secondary  flows.  Moreover,  the  geometry  is  often  asymmetric  due  to  the  volute  shape.  As 
a  result,  the  relative  motion  between  impeller  and  volute  generates  an  unstableness  which  affects  not  only 
the  overall  pump  performance,  but  is  also  responsible  for  pressure  fluctuations,  hydraulic  noises  and 
unforeseen  hydrodynamic  forces.  These  fluctuations  not  only  generate  noise  and  vibration  that  cause 
unacceptable  levels  of  stress  and  reduce  component  life  due  to  fatigue,  but  also  introduce  unfavourable 
characteristics  of  pump  performance  even  at  or  near  the  design  point.  Experimental  and  numerical 
approaches  contributed  to  the  understanding  of  the  highly  complex  flow  interactions  that  occur  in  a 
centrifugal  pump. 

Binder  et  al.  [1],  Acosta  et  al.  [2],  and  Stepanoff  [3]  conducted  some  of  the  earliest  investigations  on 
impeller  forces  in  centrifugal  pumps.  Stepanoff  proposed  a  simple  empirical  model  based  on  impeller 
geometry,  pump  head  and  capacity  to  estimate  the  radial  resultant  forces.  Agostinelli  et  al.  [4]  extended 
Stepanoff  s  model  taking  to  account  the  effect  of  specific  speed  on  radial  forces.  Biheller  [5]  developed  an 
equation  to  predict  static  radial  pump  forces  applicable  for  a  wide  range  of  pump  types  and  operating 
conditions.  Hergt  and  Krieger  [6],  Kanki  et  al.  [7]  and  Chamieh  et  al.  [8]  investigated  the  effects  of  single 
and  double  volute,  vaned  diffuser  casing  pumps  and  the  influence  of  the  impeller  blades  number  on  the 
hydrodynamic  forces.  More  recently,  de  Ojeda  et  al.  [9]  combined  the  exit  momentum  flux  and  static 
pressure  distributions  around  the  impeller  of  a  double  volute  pump  to  evaluate  a  total  resultant  radial 
thrust. 

A  number  of  authors  have  treated  the  problem  of  the  interaction  of  the  impeller  and  its  surroundings 
experimentally  (Inoue  and  Cumpsty  [10],  Sideris  and  Braembussche  [11]  and  Arndt  et  al.  [12,  13]). 

Among  others,  a  contribution  to  the  understanding  of  the  relation  between  unsteady  flow  and  mechanical 
problems  was  the  pressure  measurements  inside  a  high  specific  speed  centrifugal  impeller  operating  in  a 
double  spiral  volute  pump  using  piezoresistive  pressure  transducers  and  a  telemetry  system  (Kaupert  et  al. 
[14]). 

Instead,  Chu  et  al.  [15-17]  used  particle  image  velocimetry  complemented  with  noise  and  pressure 
measurements  to  measure  the  velocity  distribution  and  compute  the  unsteady  pressure  field  in  the  near¬ 
tongue  region  of  the  volute  of  a  centrifugal  pump.  This  pump  was  equipped  with  several  different  volutes 
having  tongue  gaps  ranging  from  7%  to  28%  of  the  impeller  radius.  The  information  available  from  PIV 
contributed  to  a  better  understanding  of  the  flow  in  centrifugal  machines  and  to  observe  and  study  the 
influence  of  the  blade  passing  in  front  of  the  volute  tongue  addressing  both  flow  structure  and  turbulence 
modelling  issues. 
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Paone  et  al.  [18]  at  VKI  used  particle  image  velocimetry  to  measure  the  flow  field  inside  the  vaneless 
diffuser  of  a  centrifugal  shrouded  pump.  They  compared  the  PIV  measurements  and  the  corresponding 
LDV  measurement  data  and  found  that  they  were  different  in  the  wake. 

In  addition  to  experimental  investigations,  numerous  analytical  studies  were  undertaken  to  predict  the 
radial  hydraulic  forces  in  centrifugal  pumps,  Domm  and  Hergt  [19],  Lorett  and  Gopalakrishnan  [20],  and 
Fongang  et  al.  [21]  to  name  only  a  few. 

Some  of  the  studies,  e.g.,  Hillewaert  and  Van  den  Braembussche  [22],  considered  the  flow  as  inviscid  and 
some  authors,  e.g.,  Longatte  and  Kueny  [23],  use  a  two-dimensional  model.  Improved  computational 
algorithms  as  well  as  hardware  development  have  contributed  to  enhance  CFD  capability.  It  is  now 
feasible  to  use  CFD  codes  for  a  realistic  prediction  of  the  complex  three-dimensional  turbulent  flow  in  the 
entire  pump  and  perform  unsteady  calculations  see,  for  example,  Zhang  et  al.  [24]  and  Gonzalez  et  al. 
[25,  26]. 

In  a  diffuser  pump,  the  centrifugal  impeller  interferes  with  its  successive  diffuser  vanes  and  produces 
pressure  fluctuations  downstream  of  the  impeller.  In  the  case  of  a  small  radial  gap  between  the  impeller 
and  diffuser  vanes,  the  magnitude  of  these  pressure  fluctuations  may  become  as  large  as  the  total  pressure 
rise  across  the  pump  (Arndt  et  al.,  [27,  28];  and  Tsukamoto  et  al.  [29]).  Dring  et  al.  [30]  indicated  two 
distinct  mechanisms  of  rotor- stator  interaction:  (I)  wake  interactions  and  (II)  potential  interactions. 
Potential  interactions  are  induced  by  inviscid  interaction  due  to  the  relative  motion  between  rotor  blades 
and  stator  vanes.  Wake  interaction  originated  from  the  impingement  and  convection  of  wakes  shed  from 
the  impeller  passages  and  moving  through  the  successive  diffuser  passages.  Some  experimental 
contributions  to  wake  interactions  in  pumps  may  be  attributed  to  the  extensive  PIV  measurements  by 
Dong  et  al.  [31],  and  Akin  et  al.  [32]  G.  Wuibaut  et  all  [33],  and  hot  wire  measurement  by  Ubaldi  et  al. 
[34]. 

However,  it  is  difficult  to  understand  the  flow  phenomena  due  to  impeller  diffuser  interaction  from  only 
experimental  studies  because  of  the  complicated  flow  structures  in  centrifugal  pumps.  Qin  and  Tsukamoto 
[34]  calculated  an  unsteady  flow  caused  by  impeller  diffuser  interaction  in  a  diffuser  pump  with  a 
singularity  method.  Following  their  works,  Shi  and  Tsukamoto  [35]  calculated  the  pressure  fluctuations 
downstream  of  the  diffuser  pump  impeller  using  two-dimensional  and  three-dimensional  unsteady  RANS 
code  with  standard  k-z  turbulence  models.  Wang  and  Tsukamoto  [36]  used  an  advanced  vortex  method  to 
calculate  the  impeller-diffuser  interaction,  in  which  the  changing  operating  points  of  pump  were  taken  into 
account. 

Unsteady  phenomena  in  diffuser  pumps  become  more  complicated  at  off-design  operating  conditions.  At 
reduced  flow  rates,  the  flow  rate  and  pressure  of  a  pump  become  increasingly  unstable.  When  substantial 
flow  fluctuations  are  propagating  at  a  low  frequency  along  the  circumference,  but  are  limited  to  a  part  of 
the  components  (e.g.,  rotor,  diffuser,  or  volute),  the  phenomenon  is  typically  referred  to  as  rotating  stall. 
Yoshida  et  al.  [37]  investigated  the  rotating  stall  instability  in  a  seven-bladed  centrifugal  impeller  with  a 
variety  of  diffusers. 

However,  the  knowledge  about  the  unsteady  pressure  fluctuations  and  the  unsteady  blade  loading  is  still 
not  satisfying.  Furthermore,  the  design  of  the  centrifugal  pumps  has  already  reached  a  level  that  only 
through  a  detailed  understanding  of  the  internal  flow  an  increase  of  the  overall  performance  can  be 
achieved.  Due  to  the  curved  passages  inside  the  impeller  and  the  volute  the  flow  is  to  be  considered  as 
three-dimensional.  Additionally,  since  the  flow  following  blade  passages  as  well  as  the  volute  casing 
interacts  with  viscous  boundary  layers,  secondary  flows  are  generated.  Therefore,  a  correct  simulation  of 
the  impeller/volute  interaction  requires  the  simultaneous  solution  of  the  three-dimensional  unsteady 
Navier-Stokes  equations  in  both  the  impeller  and  volute. 
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2.  NOMENCLATURE 
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Cp  =  (p-p2)/(pC^/2) 
Cu 
Cm 
g 

Hin 

h0 

Q 

r 

U 

Pv 

p 

c 
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area 

width 

pressure 

Absolute  Velocity 
pressure  coefficient 
absolute  tangential  velocity 
absolute  meridional  velocity 
gravity  acceleration 
input  head 

Busemann’s  head  coefficient 

flow  rate 

radius 

peripheral  velocity 
volumetric  efficiency 
density 

blockage  coefficient 
circumferential  coordinate 


Subscripts 

1  at  impeller  inlet 

2  at  impeller  outlet 


3.  VOLUTE 

In  a  turbomachine  the  impeller  does  all  the  work  while  the  volute  appears  merely  to  collect  the  flow  from 
the  runner  and  guide  it  to  the  discharge. 

It  is,  however,  wrong  to  conclude  that  the  two  parts  operate  independently  of  each  other.  The  pump 
volute  determines  the  surroundings  in  which  the  impeller  operates,  and  it  can  have  a  profound  effect  on 
impeller  performance  and  can  cause  the  impeller  to  work  inefficiently. 

The  impeller  of  a  centrifugal  pump  can  be  simply  regarded  as  a  source  of  radial  flow  on  which  is 
superimposed  a  vortex.  Therefore,  the  flow  emerges  from  the  outlet  with  a  spiral  path  whose 
circumferential  components  of  velocity  is  related  to  the  head,  while  the  meridional  component  is  related  to 
the  flow  rate.  Thus,  the  angle  of  the  spiral  flow  depends  on  the  point  of  operation  on  the  impeller 
characteristic. 


It  can  easily  be  shown  that  the  pressure  around  the  impeller  is  uniform  only  when  there  is  a  free  vortex 
velocity  distribution  in  the  volute.  For  a  perfect  fluid  the  pressure  in  all  the  points  on  a  streamline 
emerging  from  the  impeller  (fig.  1)  is  related  to  the  pressure  at  the  impeller  outlet  by  Bernoulli’s  equation: 


P  +  “pC2  —  P2  +~P^2 

and  to  the  pressure  at  B  (fig.  1)  by  the  equation  of  radial  forces: 

dp  C2 

-T  =  p — 
dr  r 


(1) 


(2) 


RTO-EN-AVT-1 43 


6-3 


Impeller  Volute  and  Diffuser  Interaction 


Figure  1:  Sketch  for  Force  Equilibrium  in  Volute  Flow. 


If  the  pressure  around  the  impeller  is  constant  these  two  equations  may  be  combined  and  lead  to  the  free 
vortex  velocity  law: 


rCu  =  r2Cu2  =  const 


(3) 


If  it  is  assumed  that  the  volute  creates  no  disturbances  to  the  flow  the  mean  velocity  in  the  throat  is 


27rr2b2C2Cm2r| 


'  Throat 


A 


Throat 


A 


1  rrMax 

Throat  12 


Throat 


bC„  dr 


r2 C  u2_  rMax 

Jro  r 


A 


Throat 


dr 


(4) 


Introducing  the  input  head  by  the  Euler’ s  formula,  the  flow  in  the  volute  can  be  related  to  the  head  by  the 
following  “volute  characteristic”: 


CU2  _  2ttb2^2r|v  Cm2 
U2  l"““2dr  U, 

Jr2  r 

gH,„  Q  1 
u|  p-  b  ^  r,U2 
Jr2  r 

It  describes  the  flow  rate  which  would  best  fit  the  volute  geometry  at  any  given  head,  and  at  a  constant 
speed,  it  represents,  in  the  H/Q  plane,  a  straight  line  through  the  origin. 


6-4 


RTO-EN-AVT-1 43 


NATO 

OTAN 


Impeller  Volute  and  Diffuser  Interaction 


The  pump  impeller  has  also  a  characteristic  relating  input  head  to  flow: 


gHin 


ho 


Q 

27tr2b2^2rivU2 


\ 

cot  |32b 

J 


(7) 


and  the  intersection  point  of  the  two  is  the  only  one  that  satisfies  the  requirement  of  both  the  impeller  and 
the  volute. 

Any  departure  from  this  condition  produces  in  the  volute  circumferential  pressure  gradients  at  the  impeller 
periphery  so  that  the  flow  in  each  blade  passage  pulsates  as  the  impeller  rotates  in  the  volute.  This  causes  a 
variation  in  the  total  head  generates,  a  subsequent  mixing  of  different  total  energy  streams  with  energy 
losses  and  flow  separation  at  the  cutwater  and  reverse  flow  from  the  volute  back  into  the  impeller  may 
take  place. 

For  evaluating  the  impeller  volute  interaction,  a  shrouded  pump  with  an  outlet  radius  of  167  mm,  an  outlet 
with  of  8  mm,  and  seven  blades  with  a  22.5  deg  outlet  lean  angle,  measured  from  the  radial  direction,  will 
be  considered.  The  single  volute  casing  is  vaneless.  The  shape  of  the  volute  casing  was  designed 
according  to  the  theory  of  a  constant  average  velocity  for  all  sections  of  the  volute,  (fig.  2).  The  gap 
between  the  tongue  and  the  volute  inner  wall  was  7  mm  and  extended  over  the  whole  volute. 


Figure  2:  Pump  Geometry. 


The  coupling  influence  on  the  reference  pump  is  illustrate  in  fig.  3,  which  shows  mean  local  streamline 
patterns  for  three  different  flow  rate. 
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(a)  53%  QDes 


(b)  100%  QDes 

Figure  3:  Streamline  Patterns. 
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When  the  spiral  path  matches  the  volute  spiral,  as  in  fig.  3b,  there  was  little  flow  disturbance  at  the 
impeller  radius  in  spite  of  the  sudden  deflection  into  the  discharge  branch. 

At  lower  and  higher  flow  rate  the  deflection  at  the  cutwater  caused  severe  energy  losses  and  was 
responsible  of  the  radial  thrust  on  pump  impeller. 

In  fig.  4  the  measured  mean  static  pressure  and  the  average  numerical  data  are  shown  for  the  volute  design 
flow  rate. 


Figure  4:  Numerical  and  Experimental  Pressure  Variation  around  Volute  100%  QDes- 

The  mean  pressure  distribution  around  the  volute  was  not  exactly  uniform  even  at  the  design  flow  rate. 
The  variation  was  limited  to  a  sudden  upward  spike  just  before  the  cutwater  and  a  corresponding 
downward  one  just  after  it.  The  sudden  jump  was  unavoidable  with  a  tongue  of  finite  thickness  and  it  was 
also  amplified,  in  this  volute,  by  the  curvature  of  the  outlet  branch. 

The  static  and  dynamic  pressure  distribution  around  the  periphery  of  the  impeller  (figs  5  and  6)  reflected 
this  dissymmetry. 


6-6 


RTO-EN-AVT-1 43 


Impeller  Volute  and  Diffuser  Interaction 


(p-p2)/(pC2A2/2)  [-] 

10  +  90 
0  +  75 

0.60 

0.45 

0.30 
0.15 
r  o.oo 
-  0.15 

-0.30 


Figure  5:  Mean  Static  Pressure  Distribution  in  the  Mid  Section  100%  Qoes- 
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Figure  6:  Mean  Static  Pressure  Distribution  in  the  Mid  Section  100%  Qoes- 


The  small  pressure  gradient  in  the  tongue  region  at  100%  QDes  did  not  influence  in  a  significant  way  the 
circumferential  distribution  of  the  outflow  at  any  section  Q(Q).  Flow  rate  and  energy  varied  close  to  the 
cutwater  only  some  points  percents,  depending  on  the  relative  position  of  the  blade  and  the  tongue  (fig.  7). 
The  velocity  fluctuation  presented  two  distinct  peaks.  The  first  was  located  at  the  interface  between  the  jet 
and  the  wake  and  the  second  trailed  behind  the  blade.  Close  to  design  conditions  the  deviation  peaks  were 
smaller  by  10  percent  and  more  evenly  distributed  around  the  volute. 
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ORGANIZATION 


Figure  7:  Instantaneous  Vector  Plot  on  Volute  Midspan  Section.  Q/QDes=100%. 


The  time  histories  of  the  pressure  coefficient  Cp  in  control  points,  located  5  mm  upstream  from  the  trailing 
edge  in  the  radial  direction  and  at  5  mm  downstream  from  the  leading  edge,  at  the  mass  flow  rate  equal  to 
the  mass  flow  rate  at  the  design  point  are  shown  in  fig.  8. 
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Figure  8:  The  Time  Histories  of  the  Pressure  Coefficient  100%  QDes- 


At  the  impeller  outlet  due  to  the  interaction  between  the  impeller  blades  and  the  tongue  of  the  volute 
casing,  amplitudes  of  pressure  fluctuation  can  be  observed.  Since  the  flow  is  incompressible,  pressure 
fluctuations  at  the  impeller  outlet  caused  by  the  interaction  between  the  blades  and  tongue  of  the  volute 
casing,  was  reflected  upstream  to  the  impeller  inlet  (fig.  8).  It  can  be  seen  that  the  amplitude  of  the 
pressure  fluctuations  within  the  impeller  passage  grows  in  magnitude  as  the  trailing  edge  of  the  blade  was 
approached.  Comparing  the  pressure  fluctuations  at  the  pressure  side  (nodal  points  1P-2P)  with  suction 
side  (nodal  points  1S-2S)  reveals  that  the  amplitude  of  the  pressure  fluctuations  on  the  pressure  side  of  the 
blade  passage  was  larger  than  on  the  suction  side.  Due  to  the  limited  interaction  between  the  impeller 
blades  and  the  tongue  the  deviation  peaks  at  inlet,  concentrated  near  the  blades  faced  with  the  tongue,  are 
negligible  (figs  8  and  9). 
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Figure  9:  Pressure  Coefficient  100%  QDes  at  Inflow. 


Figure  10  shows  the  peripheral  distribution  of  the  pressure  coefficient,  relative  to  the  mean  static  pressure 
p2  and  flow  velocity  C2  at  the  impeller  outlet,  for  several  values  of  the  flow  rate.  The  circumferential  static 
variation  was  evaluated  in  the  volute  symmetric  plane.  The  static  pressure  around  the  volute  was  nearly 
constant  for  flow  rates  close  to  the  best  efficiency  point,  but  it  exhibited  either  a  maximum  for  low 
flowrates  or  a  minimum  for  high  flow-rates  in  the  region  0  =300-330  deg.  Also,  a  sharp  pressure  decrease 
may  be  observed  close  to  the  volute  tongue  (0=10-  20  deg). 


Figure  10:  Mean  Pressure  Coefficient  Distribution  along  the  Volute  as  a  Function  of  Flow  Rate. 
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These  no  uniform  pressure  distributions  were  related  to  the  acceleration  (high  flow-rates)  or  deceleration 
(low  flow-rates)  of  the  flow  along  the  volute  for  off-design  conditions. 

Large  radial  velocities  at  high  mass  flow  resulted  in  a  high  incident  on  the  volute  tongue  and  a  subsequent 
flow  separation  in  the  exit  pipe  (fig.  11). 


0=0°  0=15° 

Figure  11 :  Instantaneous  Vector  Plot  on  Volute  Midspan  Section.  Q/QDes=135%. 


The  large  circumferential  pressure  rise  in  the  tongue  region  at  high  flow  rate  pushed  the  fluid  backward 
through  the  tongue  gab.  This  resulted  in  a  flow  separation  in  the  exit  pipe. 

If  the  instantaneous  pressure  distribution  in  the  volute  wall  was  depicted  (fig.  12),  it  can  be  observed  how 
the  wakes  advance  around  the  volute.  The  jet  and  wake  generated  by  the  blades  (fig.  12)  are  superimposed 
to  the  average  pressure  distribution  around  the  volute. 
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Figure  12:  Unsteady  Pressure  Coefficient  Distribution  along  the 
Volute  as  a  Function  of  Impeller  Position  Q/QDes=135%. 


0=-29° 


0=-22° 


0=-15° 


0=0° 


0=7° 


0=15° 


Figure  13:  Instantaneous  Dynamic  Pressure  Coefficient  ( Dp  =  (c^,  /2p)/(po2  -P2) 
Distribution  on  Volute  Midspan  Section.  Q/QDeS=135%. 
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At  the  angular  positions  opposite  to  the  tongue  (0=180  deg),  the  blade  passage  in  front  of  the  tongue  had 
little  effect,  while  near  the  tongue,  the  blade  wave  was  more  or  less  smoothed  depending  on  the  blade 
position  with  respect  to  the  tongue. 

In  general  the  pressure  fluctuations  were  comparatively  small  in  the  best  efficiency  region.  But  for  very 
low  and  high  flow-rates,  the  maximum  range,  in  the  region  close  the  tongue  (4>=10  to  50  deg)  can  reach  a 
non-dimensional  value  of  0.0  up  to  0.9  for  high  flow-rates. 

Especially  for  off-design  conditions,  the  unsteady  static  pressure  (Fig.  14)  presented  a  modulated  pattern 
along  the  volute,  with  node  and  anti-node  positions,  which  was  a  typical  effect  of  the  superposition  of 
correlated  waves.  For  QIQDes=l00%  the  pressure  fluctuations  were  observed  to  behave  similarly  along  the 
whole  volute;  they  were  synchronized  with  the  passage  of  each  blade,  and  so  the  peak  values  of  the 
fluctuations  occurred  at  different  instants  for  each  position.  These  perturbations  resulted  from  the  no 
uniform  blade-to-blade  distribution  of  the  flow  coming  out  the  impeller,  and  their  effects  were  only  local. 
For  off-design  conditions,  the  region  of  the  volute  close  to  the  tongue  presented  a  simultaneous  evolution 
of  the  pressure  fluctuations,  with  big  amplitudes  in  the  whole  region. 


Figure  14:  Instantaneous  Pressure  Coefficient  Distribution  on  Volute  Midspan  Section.  Q/QDes=135%. 

Pressure  perturbations  due  to  the  vortex  shedding  behind  the  impeller  blades  can  be  observed  also.  In  a 
small  band  around  the  vaneless  space  between  the  impeller  blades  and  volute,  complex  pressure  structures 
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can  be  observed.  These  pressure  structures  were  well  correlated  with  the  unsteady  part  of  the  relative 
velocity  at  the  exit  of  the  impeller. 

Outside  this  band,  the  pressure  field  was  smoother  and  it  can  be  observed  that  the  pressure  unsteadiness  in 
the  impeller  was  well  synchronized.  These  global  pressure  variations  were  due  to  potential  effects. 
Potential  interactions  were  induced  by  inviscid  interaction  due  to  the  relative  motion  between  rotor  blades 
and  stator. 

Indeed,  the  pressure  pulsations  around  the  impeller  outlet  and  volute  inlet  were  the  result  of  the  strong 
interaction  between  the  blades  and  the  tongue.  Comparing  the  results  for  flow-rates  above  and  below  QDes , 
the  fluctuations  in  the  region  close  to  the  tongue  were  seen  to  be  shifted  180  deg  from  each  other.  This  was 
related  to  the  shift  of  the  stagnation  point  on  the  tongue  and  the  different  characteristics  of  the 
recirculation  flow  through  the  gap.  For  low  flow-rates  the  alignment  of  the  blades  with  the  tongue 
coincides  with  a  positive  value  of  the  pressure,  whereas  it  coincided  with  a  negative  value  for  the  high 
flow-rates.  Therefore,  the  unsteady  pressure  in  the  volute  resulted  from  the  combination  of:  (/)  the  jet- 
wake  pattern  associated  to  the  continuous  blade  rotation  around  the  volute,  and  ( IT)  the  intermittent 
interaction  of  those  disturbances  and  the  volute  tongue. 

Wake  interaction  originated  from  the  impingement  and  convection  of  wakes  shed  from  the  impeller 
passages. 

With  advancement  of  the  pressure  fluctuations  were  damped,  as  it  can  be  seen  from  Figures  13  and  14,  in 
a  cross  sectional  plane  of  the  casing  at  angular  advancement  (j)  =260°,  the  pressure  fluctuations  were  small. 
In  this  cross-sectional  plane,  the  pressure  gradients  had  only  one  island  (part  I)  and  correspondingly  the 
secondary  flow  has  only  one  vortex  (part  II). 

Both  sources  can  be  initially  considered  totally  independent,  with  no  relationship  between  their  locations 
in  the  volute,  amplitudes  or  relative  phase  delay:  the  only  common  feature  was  to  radiate  harmonic  sound 
at  the  blade  passing  frequency. 

The  two  sources  can  be  reasonably  simulated  by  a  dipole,  with  amplitude  that  increased  fast  when 
diverging  from  the  best-efficiency  flow-rate.  The  positive  or  negative  combination  of  that  sound  field  with 
the  fluctuations  associated  to  the  continuous  rotation  of  the  blades  leads  to  the  modulated  pattern  in  the 
pressure  amplitude  suggested  by  Fig.  12  for  off-design  conditions. 

Figure  15  shows  the  instantaneous  pressure  field  in  control  points,  located  5  mm  upstream  from  the 
trailing  edge  in  the  radial  direction  (2P  2S)  and  at  5  mm  downstream  from  the  leading  edge  (IP  IS),  at  the 
mass  flow  rate  equal  to  135%  of  the  mass  flow  rate  at  the  design  point. 

The  potential  effect,  when  the  impeller  trailing  edge  passed  in  front  of  the  volute  tongue  was  remarkable. 

It  can  be  seen  that  the  amplitude  of  the  pressure  fluctuations  within  the  impeller  passage  grows  in 
magnitude  as  the  trailing  edge  of  the  blade  was  approached.  The  pressure  fluctuations  at  the  impeller 
outlet  caused  by  the  interaction  between  the  blades  and  tongue  of  the  volute  casing,  were  reflected 
upstream  to  the  impeller  inlet  (Figs.  15  and  16)  and  can  be  the  reason  of  a  cell  of  cavitation  located  in  the 
same  absolute  position. 
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Figure  15:  The  Time  Histories  of  the  Pressure  Coefficient  135%  Qoes- 


Figure  16:  Instantaneous  Pressure  Coefficient  Distribution  at  Inflow  Q/QDes=135%. 
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Pressure  fluctuations  at  the  impeller  inlet  and  outlet  influenced  the  mass  flow  rate  through  the  blade 
passages.  The  mass  flow  rate  through  each  blade  passage  varied  with  time  according  to  the  relative 
position  of  the  blade  passage  to  the  tongue  of  the  volute  casing  for  off-design  conditions.  The  periodic 
pressure  distribution  at  the  impeller-inlet  and  outlet  leaded  to  a  periodic  flow,  which  resulted  in  a  cyclic 
acceleration  and  deceleration  of  the  fluid  flow  inside  each  blade  passage  (Fig.  17). 


Figure  17:  Instantaneous  Deviation  from  the  Mean  Flow  Rate  and 
Total  Pressure  Elaborated  by  Each  Blade  Channel  Q/QDes=135%. 


The  diagram  shows  wide  variation  of  the  mass  and  total  pressure  values.  Near  at  best  efficiency  point  the 
circumferential  distribution  of  static  pressure  tended  to  be  uniform  and  hence  has  very  little  or  no  effect  on 
the  impeller  flow.  Therefore  at  this  operating  point,  the  circumferential  distribution  of  total  pressure,  flow 
angle,  and  total  pressure  was  unaffected  and  also  remained  uniform,  which  suggested  quite  a  rotational 
symmetric  flow  pattern. 

The  fluctuating  pressure  field  gave  rise  to  dynamic  forces  which  can  cause  fatigue  failure  of  the  pump  axis 
(especially  at  off-design  operating  conditions).  An  example  of  such  effects  can  be  seen  in  Fig.  18  where 
the  percentage  modulus  of  the  fluctuating  force  is  plotted  in  polar.  Very  large  differences  of  the  forces 
resulted  for  the  smallest  and  for  the  largest  volume  flow  of  the  pump.  At  flow  rates  in  the  proximity  of  the 
best  efficiency  point  substantially  smaller  force  differences  arose.  The  fact  that  the  curve  paths  are  not 
closed  is  related  to  the  build-up  process  of  the  transient  flow  solution  which  was  not  yet  fully  periodic. 
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Fx/F  QDes  % 

Figure  18:  Polar  Representation  of  the  Percentage  Fluctuating  Force. 


4.  DIFFUSER 

The  test  model  consisted  of  a  radial,  shrouded  impeller  with  7  blades.  The  vaned  diffuser  has  22  adjustable 
guide-diffuser  vanes  and  the  return  channel  with  1 1  continuous  vanes.  The  diffuser  allows  continuous  and 
independent  adjustment  of  the  vane  angle  and  its  relative  position  with  the  return  channel  vanes.  In  the 
present  lectures,  two  vane  angle  were  analysed  (labelled  as  Conf_l  and  Conf_2)  and  the  adjustable  guide- 
diffusers  are  faced  toward  the  return  channel  vanes. 

The  radial  gap  between  the  impeller  tip  and  the  inlet  edge  of  the  stator  vanes,  with  the  configuration  under 
test,  is  5  mm,  which  is  2.5%  of  the  impeller  radius. 

It  must  be  noticed  that  the  pump  model  has  no  volute,  but  a  collector  downstream  of  the  vaned  diffuser, 
which  allows  rather  good  axisymmetric  conditions  at  diffuser  outlet. 

Figure  19  defines  the  measuring  zone  that  was  used.  In  the  stations  SO  and  S2  mean  flow  velocity  was 
measured. 

Piezoelectric  transducers  placed  flush  with  the  wall  in  the  sections  SO  and  S2  (Fig.  19)  were  used  to  meas¬ 
ure  the  pressure  fluctuations  in  the  vaneless  gap.  Pressure  measurements  of  the  diffuser  vanes  were  made 
using  eight  equal-spaced  pressure  tapes  in  the  two  adjustable  guide-diffuser  vanes  A,  B  (fig.  20)  and  nine¬ 
teen  pressure  tapes  in  the  continuous  vanes  (blade  C  fig.  19). 
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Figure  19:  Scheme  of  the  Impeller  and  the  Vaned  Diffusers. 


Figure  20:  Pressure  Taps  Location  in  the  Adjustable  Guide-Diffuser  Vanes. 


5.  IMPELLER  DIFFUSER  INTERACTION 

In  Figures  21  and  22  the  pump  head  and  the  efficiency  are  reported  together  with  the  shaft  power  for  the 
two  adjustable  guide-diffuser  vane  angle  analysed  in  parallel  with  the  vaneless  diffuser. 
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Figure  21 :  Overall  Performance  Conf_1  and  Conf_2. 


Figure  22:  Power  Spectrum  of  Pressure  Measured  in  the  Inlet  Duct  (Sez.  SO  -  Conf_1). 


The  characteristic  of  the  pump  becomes  slightly  unstable  between  Q/QDes  =  0.55  to  0.70  for  both  ad¬ 
justable  guide-diffuser  configurations. 

Comparing  the  curves  obtained  with  the  two  configurations,  the  reduction  of  the  guide-diffuser  vanes 
angle  moderated  the  intensity  of  the  instability  but  does  not  remove  it. 

Furthermore  the  maximum  efficiency  is  decreased  slightly  and  it  is  moved  towards  lower  flow  rate. 

Figs.  22,  23,  24  and  25  show  the  pressure  powerspectra,  measured  for  the  two  diffuser  configuration. 
Pump  stator-rotor  interaction  mechanisms  generated  both  discrete  frequency  noise  and  broadband  noise. 
The  discrete  frequency  components  of  the  pressure  fluctuations  in  the  impeller  consisted  principally  of 
fstator=  nBiade  Diffuser  03/271  (St=l)  and  its  higher  harmonics  and  the  impeller  passing  frequency  fstator=  nBiade 
impeller  03/271  and  its  higher  harmonics  basically. 
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Figure  23:  Power  Spectrum  of  Pressure  Measured  in  the  Inlet  Diffuser  (Sez.  S2  -Conf_1). 
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Figure  24:  Power  Spectrum  of  Pressure  Measured  in  the  Inlet  Duct  (Sez.  SO  -  Conf_2). 

003i 


Figure  25:  Power  Spectrum  of  Pressure  Measured  in  the  Inlet  Diffuser  (Sez.  S2  -  Conf_2). 

The  broadband  component  was  associated  with  the  turbulent  flow.  Noise  was  generated  as  the  inertia  of 
vortical  disturbances  was  modified.  This  occurred  essentially  due  to  the  sudden  change  in  the  boundary 
conditions  applied  to  the  vorticity  as  vortices  interact  with  solid  boundary. 
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At  high  flow  rates  the  spectral  shape  was  very  similar  for  each  discharge  configuration  suggesting  that  the 
dominant  type  of  pressure  mechanism  did  not  change  significantly  with  modifications  to  the  impeller  dis¬ 
charge  configuration. 

The  influence  of  the  interactions  on  the  pressure  becomes  much  stronger  near  the  impeller  trailing  edge. 
Fig.  26  and  27  show  the  power  spectrum  of  pressure  at  the  discharge  of  two  consecutive  impeller  blades  in 
the  two  configurations.  The  spectral  shapes  were  very  similar  for  each  discharge  configuration  and  they 
were  dominated  by  tones  of  the  diffuser  blade  passage  frequency.  The  influence  on  the  pressure  side  was 
much  bigger  than  that  on  the  suction  side.  In  that  case  the  influenced  region  on  the  suction  side  reached  up 
to  the  streamwise  location  of  about  70%  far  from  the  leading  edge,  but  that  on  the  pressure  side  can  reach 
up  to  40%.  The  influenced  regions  of  the  diffuser  and  the  return  channel  were  the  whole  lengths  of  the 
vanes. 


Figure  26:  Power  Spectrum  of  Pressure  at  the  Discharge  of  Two  Consecutive 
Impeller  Blades  (Pressure  Side  -  Conf_1  -  Q/QDes=0.62). 


Figure  27:  Power  Spectrum  of  Pressure  at  the  Discharge  of  Two  Consecutive 
Impeller  Blades  (Pressure  Side  -  Conf_2  -  Q/QDes=0.62). 
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Figure  28:  Pressure  Amplitudes  in  the  Midspan  of  the  Adjustable 
Guide-Diffuser  Vanes  A,  Conf_1,  Q/QDes=0.75. 


When,  an  impeller  blade  began  to  approach  a  diffuser  vane,  the  interactions  decreased  the  pressure  in  the 
whole  diffuser  region  (Fig.  16),  resulting  in  pressure  decrease  in  the  return  channel.  When  the  approaching 
process  finished,  the  pressure  on  the  suction  side  of  the  diffuser  firstly  decreased  till  the  streamwise 
location  of  0.75,  then  increased;  however  the  pressure  on  the  pressure  side  increased  within  the  whole 
streamwise  length.  When,  the  impeller  trailing  edge  moved  away  from  the  diffuser  leading  edge,  the 
influence  decreased  in  all  the  regions. 

The  jet-wake  discharged  from  the  impeller  sweeps  periodically  over  the  diffuser  and  induced  a  very 
strongly  varying  periodic  flow  in  the  entry  zone.  The  jet- wake  flow  emerging  from  the  impeller  gave  rise 
to  substantial  circumferential  and  spanwise  distortions  that  did  not  mix  out  before  arriving  at  the  diffuser 
vanes  entry  zone. 

At  the  nominal  operating  point,  all  the  measured  and  calculated  velocity  vectors  were  directed  from  the 
impeller  into  the  diffuser,  there  was  no  evidence  of  back  flow  and  the  pressure  fluctuations  were 
dominated  by  the  discrete  frequency  components.  At  Q=0.86QDes  some  flow  separations  were  initiated 
within  the  two  diffuser  close  to  the  hub  surface  (Figs.  29  and  30).  The  flow  pattern  within  the  continuous 
vanes  is  very  complex,  because  of  the  inlet  shear  flow  and  the  large  flow  passage  curvature  both  on  the 
meridional  and  blade  to  blade  passage  plane.  Close  to  the  hub  the  flow  follows  the  suction  surface  of  the 
diffuser  blade  at  the  inlet  but  the  flow  around  the  midpitch  location  started  reversing  upstream.  Just  before 
the  blade  to  blade  bend  a  vortex  in  the  counter  clockwise  direction  was  observed  due  to  the  overloading  of 
the  blade  (fig.  29).  In  the  adjustable  guide-diffuser  vanes  the  vortex  was  considered  to  be  created  by  the 
blockage  effect  of  the  contiguous  continuous  diffuser  systems  (fig.  29).  The  reverse  flow  was  not  axi- 
symmetric  and  was  more  extensive  between  the  midpitch  and  the  pressure  side  of  the  diffuser  blade. 
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Figure  29:  Instantaneous  Velocity 
Vectors  Close  to  the  Hub  at 


Figure  30:  Instantaneous  Velocity 
Vectors  Close  to  the  Hub  at 


Q-0.86Qoes  Conf  1 . 


Q_ 0.6/Qpes  Conf  1 . 


The  onset  of  the  recirculation  was  confirmed  by  an  increase  of  the  pressure  level  measurements  (Figs.  23 
and  25)  and  by  the  change  of  the  power  spectrum  of  pressure  close  to  the  hub  at  the  discharge  of  the 
impeller  (Figs.  26  and  27)  measurements  in  the  range  between  0.00  and  0.25  kHz. 

With  diminished  flow  rates,  the  backflow  area  expanded  and  the  onset  of  the  recirculation  into  the  impeller 
from  the  diffuser  was  detected  near  the  leading  edge  of  adjustable  guide-diffuser  vanes  (Figs.  30  and  32). 


Figure  31:  Instantaneous  Velocity 
Vectors  Close  to  the  Hub  at 
Q=0.86QDes  Conf  2. 


Figure  32:  Instantaneous  Velocity 
Vectors  Close  to  the  Hub  at 
Q=0.67QDes  Conf  2. 


Two  flow  phenomena  can  be  observed  in  the  snapshot  of  the  flow  field.  First  the  onset  of  the  recirculation 
into  the  impeller  from  the  diffuser  can  be  detected  near  the  leading  edge  of  adjustable  guide-diffuser  vanes 
(Figs  30  and  32).  Second  the  start  up  of  a  clockwise  vortex  was  observed  close  to  the  shroud  at 
Q=0.67QDes-  The  development  of  the  recirculation  seems  to  be  confirmed  by  the  progressive  increase  of 
the  mean  pressure  level  noise  measured  experimentally  while  the  onset  of  a  second  reverse  area  was 
shown  by  the  change  of  the  slope  of  the  mean  pressure  level  in  the  range  between  0.00  and  0.25  kHz. 
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At  flow  rate  lower  than  Q=0.53QDes  the  recirculation  zones  close  to  the  hub  disappear  and  the  back  flow 
area  was  confined  close  to  the  shroud.  The  change  of  the  recirculation  zone  from  the  shroud  to  the  hub, 
considered  to  be  very  important  from  the  stability  of  the  head  capacity  characteristic  at  part  load  (Figs.  21 
and  22),  was  confirmed  by  the  abrupt  decrease  of  the  noise  due  to  the  disappear  of  one  of  the  noise  source 
in  the  range  between  0.0  and  3.0  St  (Figs.  24  and  25). 

Unsteady  pressure  on  the  blade  surfaces,  produced  by  impeller  blade  interactions  with  the  instability  at  the 
impeller  discharge  were  believed  to  be  the  major  source  of  noise.  In  the  noise  spectrum,  broad  humps 
appeared  which  were  centred  to  the  peaks  in  source  strength  spectral  distribution  initially  identified  at 
higher  flow  rate.  All  modes  in  the  pressure  spectrum  generated  noise.  The  peak  at  St=0.286  generated 
more  acoustic  noise  than  the  modes  with  higher  amplitudes  (Figs  24  and  25)  and  was  identified  as  a 
pressure  instability  that  formed  a  rotating  pattern  around  the  impeller. 


The  midspan  section  of  the  diffuser  vanes  seems  to  be  interested  only  marginally  from  the  recirculation 
and  the  relieves  evidenced  regular  recoveries  of  pressure  both  in  the  adjustable  guide-diffuser  and  in  the 
continuous  vanes  (figs.  33,  34,  35  and  36). 


L/Ljviax  L/Ljyiax 


Figure  33:  Average  Pressure  in  Midspan 
of  Adjustable  Guide-Diffuser  Conf  1. 
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Figure  35:  Average  Pressure  in  Midspan 
of  Continuous  Vanes  Conf  1. 


Figure  34:  Average  Pressure  in  Midspan 
of  Continuous  Vanes  Conf  2. 


Figure  36:  Average  Pressure  in  Midspan 
of  Continuous  Vanes  Conf  2. 


The  inlet  flow  spectra  showed  low  amplitude  at  the  blade  passage  frequency  and  higher  harmonics  for  quite 
all  flow  rates.  The  only  tones  found  in  the  spectra  was  related  to  the  sub-tonal  harmonics  (0<St<l)  (figs  22 
and  24). 
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